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Abstract 

Tensor models can be interpreted as theory of dynamical fuzzy spaces. In this 
paper, I study numerically the fluctuation spectra around a Gaussian classical 
solution of a tensor model, which represents a fuzzy flat space in arbitrary dimen- 
sions. It is found that the momentum distribution of the low-lying low-momentum 
spectra is in agreement with that of the metric tensor modulo the general coor- 
dinate transformation in the general relativity at least in the dimensions studied 
numerically, i.e. one to four dimensions. This result suggests that the effective 
field theory around the solution is described in a similar manner as the general 
relativity. 
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1 Introduction 



There have been various discussions concerning the fuzziness of spacetime PQ-jE]- In spite of the 
absence of fundamental theory of spacetime, it is generally argued that quantum gravitational 
fluctuations invalidate the reality of the classical notion of spacetime, i.e. a smooth continuous 
manifold in the general relativity. Therefore, it is widely believed that this classical notion 
should be replaced in some way with a new quantum notion. Another unsatisfactory feature 
of this classical notion is that a spacetime and the gravitational degrees of freedom, i.e. a 
metric tensor, are independent distinct objects, i.e. in principle, it is allowed to consider a 
spacetime without a metric tensor. This kind of insufficiency seems to exist also at least in the 
present formulation of string theory, i.e. a string and a spacetime are quite different physical 
objects, although they are tightly related in several ways, such as the restrictions of spacetime 
dimensions. Therefore it is desirable that a new notion of spacetime contains a spacetime itself 
and the physical degrees of freedom in an inseparable fashion. The construction of a quantum 
theory of spacetime or quantum gravity would essentially be equivalent to finding such a new 
notion of spacetime which contains inseparably the general relativity in its classical limit. 

An interesting candidate for such a new notion of spacetime comes from fuzzy spaces. The 
central idea is the equivalence between a spacetime and the algebra of functions on it. For a 
usual space, this is the implication of Gelfand-Naimark theorem [9j . Considering a deformation 
of an algebra of functions on a usual spacetime, a fuzzy spacetime can be obtained. When 
the deformed algebra is noncommutative, it defines a noncommutative spacetime [9]-[12j. 
The deformed algebra can even be nonassociative [13]- [2D]- This direction of deformation 
to nonassociativity would be important in allowing an algebra to contain the information of 
geometric structures of a fuzzy spaceB 

In view of the above idea of fuzzy spaces, a dynamical theory of spacetime may be obtained 
by a dynamical theory of fuzzy spaces. A fuzzy space is defined by an algebra, and an algebra 
can be given by a rank-three tensor C a b C , which defines the multiplication rule among the 
algebraic elements f a as f a fb = C a bf c - Therefore it would be possible that a dynamical 
theory of a rank-three tensor C a b C , i-e. a tensor model, can be regarded as a dynamical 
theory of spacetime. This direction was pursued by the present author [2T]-[23], and a number 
of physically interesting classical solutions of tensor models were found numerically, which 
represent commutative nonassociative fuzzy spheres and tori in various dimensions [23J. 

It is interesting to note that tensor models have also appeared in some other contexts of 
quantum gravity. Tensor models were considered as analytic expressions describing the dy- 



*This was partially insisted in the previous work [23j by the present author. In noncommutative spacetime, 
the geometric properties are contained in the Laplacian rather than in the algebra [9]. 
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namical triangulation model of spacetime in more than two dimensions and topological lattice 
gravities [24j-[28j. They also appeared in relation with the loop quantum gravity [29]-[3T]. 
The main difference of these previous approaches from the present one is the interpretation of 
tensor models. In the previous approaches, the Feynman graphs of the amplitudes of tensor 
models correspond one-to-one to the dual diagrams of triangulated manifolds. Computations 
of tensor models under this interpretation need essentially quantum treatment of tensor mod- 
els, which, however, requires future technical developments. On the other hand, under the new 
interpretation, a classical solution corresponds to a background spacetime, and the fluctuations 
around it will correspond to the fluctuations of the physical fields on the spacetime. 

In this paper, I will numerically study the fluctuation spectra around a classical solution of 
a tensor model. This solution has a simple Gaussian form, and represents a commutative but 
nonassociative fuzzy flat space in arbitrary dimensions [T5]. The main result of this paper is 
that the number of the low-lying low-momentum fluctuation spectra in each momentum sector 
agrees exactly with that of a metric tensor modulo the general coordinate transformation in 
the general relativity. This suggests the possibility that the effective field theory around 
the classical solution can be described in a similar manner as the general relativity. If this 
is definitely proved, this specific tensor model contains both the spacetime and the metric 
degrees of freedom in an inseparable fashion, and can be qualified as a quantum gravity. 

This paper is organized as follows. In the following section, a tensor model is defined, which 
has a one-parameter family of a Gaussian classical solution. In Section [3j the explicit procedure 
to obtain the fluctuation spectra and the zero modes generated from a symmetry breaking by a 
classical solution are discussed. In Section HI discretization of momenta and momentum cutoff 
are introduced as regularization for numerical study, and numerical solutions very similar to 
the Gaussian solution are obtained in one and two dimensions. The numerical analysis of the 
fluctuation spectra around the numerical solutions is carried out. It is shown that there exist 
some non-zero but very light modes, the distribution of which has a characteristic feature 
depending on the dimensions. In Section [5j the distribution of the very light modes is studied 
further by an approximate but more efficient numerical method in one to four dimensions. In 
Section El the characteristic feature of the distribution of the very light modes is shown to be 
explained in a simple manner, if they are identified as the metric in the general relativity. The 
final section is devoted to the discussions and future prospects. 

2 A tensor model with a Gaussian solution 

The tensor model proposed in [21] has a real three-tensor C a t, c and a real symmetric two- 
tensor g ab as its dynamical variables. The model has an invariance under a general linear 
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transformation, which can naturally be identified as a fuzzy analog of the general coordinate 
invariance in general relativity under the interpretation of the tensor model as a dynamical 
theory of fuzzy spaces. Without any simplifications, however, the model is quite complicated 
in general. Therefore, in this paper, it is assumed that the two-tensor g ab be non-dynamical, 
and also that the three-tensor C a b c be symmetric under the permutations of the indices. The 
former assumption breaks partially the general linear symmetry, and the model given below 
has only an orthogonal group symmetry. Thus the assumption seems to break the above 
connection of symmetry between the tensor model and the general relativity, but the result 
of this paper will suggest that they are actually intimately related. The reason is not clear 
presently, but it may be that the general relativity is related to a broader class of tensor 
models than expected solely from symmetry. 

The purpose of this section is to construct a tensor model which has a classical solution with 
a simple Gaussian form. This classical solution represents a fuzzy space which is essentially 
the same one introduced in [19J to study a field theory on a nonassociative spacetime. In the 
momentum basis, the fuzzy space is defined by the algebra f pl f P2 = Cp lP2 P3, f PZ with 

C PlP2P3 = A exp (-a (p\ + p\ + pf) ) 5 D (pi + p 2 + p 3 ) , 

gPV* = 5 D ( Pl + P2 ), (1) 

where A and a are positive constants, D is the dimension of the space, and PiS are continuous 
D dimensional momentum vectors. Only the Euclidean signature is considered in this paper. 
Here I have used the notation C a ^ c with a bar for the Gaussian solution to distinguish it from 
the dynamical variable C a bc- Similar notations will be used for some other tensors which will 
be defined below. Note that g ab is non-dynamical and takes the fixed value ([1]). In a coordinate 
basis, after Fourier transformation, ([1]) becomes 

C Xl x 2 x 3 ~ exp (-/3 ((xi - x 2 f + (x 2 - x 3 ) 2 + (x 3 - xi) 2 )) , 
g>™ ~ 5 D (x 1 -x 2 ), (2) 

where f3 ~ I /a. The reality condition on the dynamical variable C a bc is imposed in this 
coordinate representation. 

The multiplication algebra f Xl f X2 — C XlX2 X3 f X3 defined by (T5]) is commutative but nonas- 
sociative. This fuzzy space is a fuzzy analogue of a D dimensional flat space. In fact, (DO) 
(or ([2])) has obviously the Poincare symmetry, and, in the a — > {(3 — > oo) limit with fixed 
coordinates^, the multiplication algebra approaches f Xl f X2 ~ 5 D (xi — x 2 ) f X2 , which is the 
algebra of the functions on a usual continuum flat space by the identification f Xl = 5 (x — x±). 

^Whilc the all-over scale has no meaning, since it can be rescaled by the rescaling of a, the relative scale 
has a definite meaning. 
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Figure 1: From the left to right, K a b, Ji a bc and I a b c . The three-vertex represents the three- 
tensor C a b c and the connecting lines represent the contractions with g ab . The open ends are 
associated with the indices of the tensors. 

The nice property of the Gaussian form ([I]) is that the contraction of the indices can 
explicitly be computed, and the resultant tensor has again a Gaussian form. The integration 
measure of momentum in a contraction is defined by dPp. One can easily compute the 
following quantities shown graphically in Figure [TJ 
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exp (-5a (pI+pI+pI)) S D {pi +p 2 +P3)- (5) 



Let us find an equation of motion satisfied by (Q. From the observation of the decay 
exponents in (j3J) and (jSJ), one finds that the quantity defined by 

W PlP2P3 = C PlP2P3 - K p ^K p ^K p * pl J PlP2P3 , (6) 

where K~ b x is the inverse of K ab , vanishes for (CQ), if A satisfies 
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Thus 

W P1P2P3 = (8) 

is an equation of motion which has (pQ) as its solution. Note that (jHJ) has the solution (flj) 
irrespective of D, and gives a kind of unified description. 

There exist various actions which derive the equation of motion (jSJ). The simplest one 
would be 

S = g aa 'g bb 'g cc 'w abc W a , hlc ,. (9) 

Since the action is positive definite and the action vanishes at the solution (pQ), the solution is 
stable, i.e. all the fluctuation spectra are non-negative. 



3 Fluctuation spectra around a classical solution 

In this section, I will discuss the explict procedure to obtain the quadratic fluctuation spectra 
of the tensor model around a classical solution, and also the zero modes generated from a 
symmetry breaking. 



3.1 The quadratic fluctuation spectra 



The second derivative of S with respect to C abc at a solution C abc = C® bc to W abc = is given 
by 

g g g JJ Br Br — 

C=C° UK ^abc V^a'b'c' 



d 2 S 



(10) 

c=c° 



2 dC abc dC a > b > c > 

However, this expression is not appropriate for obtaining the fluctuation spectra. This is 
because the derivative should be taken with respect to the independent normalized variables, 
but not all of the components of C abc are independent due to the permutation symmetry of its 
indices. In the followings, let me assume that the range of an index of the tensor is finite, i.e. 
a = 1, 2, • • • , n. This is not true for the continuum momentum index in the preceding section, 
but a regularization will be introduced later for the numerical computations. Then a natural 
normalization can be defined by the following SO(n,R) invariant measure, 

(5C) 2 = g aa 'g bb 'g cc '5C abc 5C a , b/c , 

= rf (Pl'P2,P3) l<5C(pi, P2 ,P3)| 2 - (11) 

(Pl,P2,P3) 

Here the summation is over the momentum set (j)i,P2,P3), where the order of the momenta 
Pi)P2,P3 is not cared, i.e. (pi,P2,P3) = (P2,Pi,P3) — 4 • • , and d(pi,p2,P3) denotes the number 
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of times the set (pi,P2,P 3 ) appears in the simple summation ^2 pip2P3 - Namely. 
d(p l7 p 2 ,p 3 ) = 

Therefore the normalized independent variables can be taken as 



1 for pi = p 2 = p 3 , 

3 for pi = p 2 ^ p 3 , P2=P3l L Pi, or p 3 = p 1 ^ p 2 , 
6 otherwise. 



a 
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Thus the fluctuation spectra should be evaluated by 
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where the partial derivatives in the last line are taken as if all the C a bc are independent, and 
the sum over a, a' is over all the permutations of 1,2,3 to compensate this virtual treatment. 
The fluctuation spectra can be obtained by diagonalizing the matrix (TBI) . 

Since W a bc is a complicated function of C a b c , the computation of the explicit formula of the 
partial derivatives in (Tl"4"l) is straightforward but lengthy. All the necessary formula are shown 
in Appendix A. 



3.2 The zero modes from the symmetry breaking SO(n, R) / SO(2) 



D 



In this subsection, I will study the zero modes of the fluctuation spectra. A number of zero 
modes appear, since a classical solution breaks the SO(n, R) symmetry of the tensor model in 
general. Let me start with the standard discussion. 

Suppose there is a classical solution C a b c = C® bc , 

U^abc C=C° 

The SO(n, R) symmetry of the action implies that, a new solution can be generated by apply- 
ing a broken transformation to C® bc . By considering an infinitesimally small transformation 
of (|T5l) . this implies 

(M/C° d , ef + M/Cl f + M/C° defl ) = 0, 



dC a bcdCdef 



(16) 



c=c° 
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where M a a is a broken generator. Namely, 

8C& = M/C° a , bc + M b b 'C° ab , c + M/C° abc , (17) 

is a zero mode of the quadratic fluctuation spectra. Thus, in general, there is a one-to-one 
correspondence between the broken generators and the zero modes of the fluctuation spectra 
around a classical solution. 

The general coordinate transformation of the metric tensor, 5g fiu = V^tv + V^t^, in the 
general relativity is not linear in g^ v . There is an old idea that the non-linear transformation 
comes from a spontaneous breaking of the local translational symmetry [32j, in analogy with 
the standard non-linear realizations of global symmetries. The present setting of the tensor 
model fits very well with this old idea, since the symmetry of the tensor model can naturally 
be identified with the fuzzy analogue of the general coordinate invariance [21] and a classical 
solution breaks the symmetry in general. Then the zero modes of the fluctuation spectra 
should be identified with these gauge zero modes, which are unphysical. 

In the numerical investigations in the following sections, the values of momenta will be 
discretized in integers, and a cutoff will be introduced. Namely, an index of the tensor takes 
the form, 

a = (p\p 2 ,--- ,p D ), (p^-^-L + l,--- ,L). (18) 

I will also make an ansatz that there remains a D-dimensional translational invariance SO(2) D 
on a solution. These assumptions mean that the background space corresponding to a solution 
is actually a fuzzy D-dimensional flat torus. The broken symmetry is given by SO(n, R) / SO(2) D , 
and thus the total number of the zero modes is given by 

TliTl — 1 ) 

#zero totai = V 2 1 - D, (19) 

where n is the size of the range of the index ( fl8|) . i.e. n = (2L + 1) D . 

The above counting of the zero modes can be made more detailed by counting them in 
each momentum sector, which will be used in numerical investigations in the following sections. 
The generator of SO{2) D is given by 

M; iP2 = i(v- Pl )5 pl _ P2 , (20) 

where v is a real D-dimensional vector. It acts on C a b c as 

iv ■ (p! + p 2 + ps) C P1P2P3 . (21) 

Therefore, because of the translational symmetry, the tensor of a solution C p p3 takes non- 
vanishing values only at the vanishing momentum, pi + P2 + P3 = 0. This implies that 
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the momentum of a zero mode SC^ C in (|17p comes solely from that of a generator M a &. 
Thus counting the zero modes in each momentum sector is equivalent to counting the broken 
generators in each momentum sector. 

There are two conditions for M ab to be a generator of SO(n, R) as 

M P1P2 = -M P2P1 , (22) 

where * denotes the complex conjugate. The first condition comes from the reality of the 
generator in the coordinate representation^, and the second one comes form the anti-symmetric 
property of the indices of the generators of SO(n,R). The total momentum of M piP2 is 

P = Pl +P2- 

When the total momentum vanishes, i.e. p — 0, the conditions (T2"2"l) are given by 

M_ P1P1 = -M P1 _ P1 , (23) 

where p\ is an integer vector as ffl8|) . These conditions imply that M^^ must vanish, and for 
Pi 7^ 0, M pi _ pi are pure imaginary and only half of the generators are independent. Therefore 
the number of the independent generators in the p = sector is given by 

#M(3) = i(n-l). (24) 

When the total momentum p does not vanish, the conditions become 

Ml pi+ppi = M pi _p_ Pl , (25) 
M- Pl+P Pl = —M Pl - Pl + P , (26) 

where both p\ and —p\ + p must be integer vectors as ffl8l) . The second condition implies 
that, when there exists a p\ which satisfies 2pi = p, i.e. when p is an even vector, M PlPl must 
vanish. It also implies that, for 2p% ^ p, only half of M_ P1+P Pl are independent, each of which 
is a complex number. The first condition implies that the generators in the momentum — p 
sector is determined by those in the momentum p sector. In fact, under the reality condition 
in the coordinate basis, the modes must be a combination of those in the p and — p sectors, 
in the same was as sin(px) = (e tpx — e~ ipx )/(2i) or cos(px) = (e ipx + e~ vpx )/2. Therefore this 
effectively generates an additional factor 1/2 to each p and — p sectoJ!. Counting the number 



■^Note that the reality condition is imposed in the coordinate representation like @. 
§This rough treatment can be justified in the actual usage in the following numerical study. 



8 



of the possible p\ under the condition that both pi and —pi + p must be in the range of f[TB"j) . 
and subtracting 1 when p is an even vector, one finally obtains the number of the independent 
generators of SO(n, R) in the momentum p sector as 



#M(p) = M J] (2L + 1 - bi) - even(p) 

\i=l / 



(27) 



where is the z-th component of the momentum p, and 



evenf ) — I ^ ^ or ^ even vector, 
1 otherwise. 

Note that the result (j2"7|) includes also the special case (1241) . Note also that, since the momen- 
tum of an index is bounded as (118]) . each component p % of the total momentum p = p\ +P2 in 
(127]) is bounded by < 2L. 



The unbroken genenerators (1201) are in the p = sector. Subtracting this from ( 1271) . one 
obtains the formula for the number of the zero modes in each momentum sector as 

# Z ero(p) = { i (Hi. < 2i + 1 - \P'\) " —(f)) " ™* ta \P'\ * 2i - (29) 
otherwise. 



4 Fluctuation spectra around the Gaussian-like numer- 
ical solutions in one and two dimensions 



The Gaussian solution ([T]) is a function of continuous infinite momenta, and is not suited for 
direct numerical investigations. Therefore (Tl8l) is introduced as regularization. Then, however, 
(CQ) does not satisfy the equation of motion (jSJ). There are two ways to deal with this problem. 
One is to obtain numerical solutions similar to (p]) to the equation of motion ([8}, and study 
numerically the fluctuation spectra around it. The other is to use ([1]) as an approximate 
solution to the equation of motion (jSJ), and evaluate numerically the spectra by substituting 
(PQ) directly into ([141) . The former way will be carried out in this section, and the latter will 
be the subject in the next section. The former way is definite but takes longer computational 
time, and will be carried out only in dimensions D=l,2 with smaller L. The latter way is an 
approximate treatment but is more efficient, since it does not contain the process of obtaining 
numerical solutions, and will be carried out in D = 1, 2, 3, 4 for larger L. 

The numerical computation is carried out on a Windows XP 64-bit workstation with two 
AMD Opteron 275 processors (dual core, 2.2GHz). The total memory size is 8 GB. The 
main numerical computation is carried out with Intel C++ compiler 9.1 and 10.0 with Open 
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Figure 2: The values of for D = 1, L = 3, 5, 10 (left to right). The horizontal axis is p. 

MP parallelization. AMD core math library (ACML) 3.6.0 is used for some basic routines. 
Some supportive computations such as drawing graphs and cross checks are carried out with 
Mathematica 5.2 personal grid edition. 

Let me first consider the D = 1 case. Let me make some ansatz in obtaining the solutions to 
((El). One is the translational invariance of the solution, namely, C pip2P3 ^ only at P1+P2+P3 = 
0. This was assumed in counting the zero modes in Section [3721 Another is that C pip2P3 be 
real. The last one is 

C = f (30") 

VlVlPZ ~Pl-P2-P3> \ J 

which is a reflection symmetry of the corresponding fuzzy S 1 . 

The numerical solutions are obtained by taking starting values of C PlP2P3 as (pQ) with a ~ 
1/L 2 , searching for the minimum of the sum of the equation of motion (JS]), 

f(C)= £ W* 1P2P3} (31) 
(P1P2P3) 

and checking whether the minimum can reasonably be regarded as vanishing within the com- 
putational accuracy. The Nelder-Mead method is used for the optimization. 

The numerical computation is carried out for L — 3,5, 10. The raw numerical results are 
shown in Appendix [B] In Figure [2j the values of K v p , which is defined by neglecting the bar 
in ([3]) (or in Figured]), are plotted for the numerical solutions. The shapes look very much 
like Gaussian, and therefore the numerical solutions can be regarded as the correct discrete 
solutions to ([8]) corresponding to the continuum Gaussian solution ([1]). 

Let me first look at the L = 3 case. The spectra in Table [2] contain values < 10 -14 . The 
number of these tiny spectra in each momentum sector agrees with that of the zero modes 
from the formula ( 1291) . Therefore these modes can be identified with the zero modes generated 
by the symmetry breaking SO(n, R)/ SO (2). Another thing one notices is that there exists a 
non-zero but very light mode with ~ 2 x 10~ 3 only at the p = sector, while the others are 
> O(0.1). The qualitative features of the spectra for L = 5, 10 are similar with L = 3, i.e. the 
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Figure 3: The fluctuation spectra for D — 1, L — 3, 5, 10 (left to right). The horizontal axis 
is the momentum p, and the vertical one is the spectrum. 

agreement with (129]) . the existence of a very light mode only at the p = sector, and that the 
others are > O(0. 1). The very light modes have the spectra ~ 1 x 10~ 6 and ~ 3 x 10~ 7 for 
L = 5, 10, respectively. The spectra are plotted in Figure EJ 

The D = 2 case is carried out essentially in a similar manner. To reduce the number of 
the variables in searching for a solution, some reflection symmetries are imposed. One is the 
reflection symmetry in each direction, 

(32) 

The other is a symmetry which exchanges the two coordinates, 

C (p\,pI) {pIpD (p\,pD = C °{pIv\) (pIpD (pIpIY ( 33 ) 

The meaning of these symmetries is that the space is a fuzzy analogue of S 1 x S l , where the 
two S l, s are congruent and perpendicular to each other, and each of them has a reflection 
symmetry. 

The numerical computation is carried out for L = 2,3. Numerical solutions similar to 
([1]) are found. In fact, K^'^ for the numerical solution for L = 3 is plotted in Figure HI 
the profile of which is similar to Gaussian. The raw numerical results are summarized in 
Appendix [B] The spectra at low momentum sectors for L = 3 are plotted in Figure [5j Some 
qualitative features are similar to the D = 1 case. The tiny spectra can be identified with 
the zero modes, the number of which agree with fl29|) . There exist also non-zero but very 
light modes, the distribution of which, however, is different from the D = 1 case. There exist 
three very light modes at the p = (0, 0) sector, and one very light mode at each non-vanishing 
momentum sector. For L = 3, the spectra of these very light modes are 1.3 x 10~ 12 , 3.0 x 10~ 4 
and 1.1 x 10~ 3 at the p = (0,0) sector, and 2.1 x 10~ 3 , 9.6 x 10~ 4 and 6.9 x 10" 3 in the 
p = (1, 0), (1, 1) and (2, 0), respectively. These modes are clearly distinguished from the other 
heavy modes with spectra > O(0. 1). In the larger momentum sectors, however, the hierarchy 
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Figure 4: The profile of K^^J for D = 2, L = 3. The horizontal axes are (p 1 ,^ 2 ). 
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Figure 5: The fluctuation spectra for D = 2, L = 3. 
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between the very light modes and the other heavier modes cannot be clearly observed. This 
situation is expected to be improved in larger L. 

The number of the dynamical variables C a b c expands very rapidly with the increase of 
D, L, and it seems very hard to carry out the same computations for larger D, L in the present 
computational facilities. Our final goal is of course D = 4. Therefore, in the next subsection, 
approximate treatment will be carried out for larger D, L including D = 4. The numerically 
correct results in this section will be used to support the results of the approximate treatment 
in the following section. 

5 Distribution of the very light modes with an approx- 
imate method 

The results in the previous section show that the spectra can be classified into three cat- 
egories. One consists of the heavy modes with spectra > O(0. 1), another consists of the 
very light modes, and the other consists of the zero modes from the symmetry breaking 
SO(n, R)/SO(2) D . The distribution of the very light modes seems to be characteristic of the 
dimensions D. The main purpose of the present section is to study numerically the distribution 
of the very light modes in one to four dimensions by using a more efficient but approximate 
method. As has been checked for some D, L in the preceding section, I assume that the number 
of vanishing spectra is given by (129]) also for other D, L. Namely, zero modes are in one-to-one 
correspondence to the broken generators of the symmetry breaking SO(n,R)/SO(2) D . 

An inefficient part of the numerical computations in the previous section is the process to 
find the numerical solutions by minimizing /(C). The reasons for the inefficiency are not only 
the complexity of the equation of motion ([H]), but also the existence of the very small spectra 
at the p = sector. The function f(C) has a very gentle slope in this direction, and it is hard 
to find its minimum. This time-consuming step can be skipped, if the analytical expression 
([T]) can be used as an approximation to a solution to the equation of motion The errors 
of the approximation come from the cutoff and the discretization of momentum introduced in 
fTTS]) . One may appropriately choose the parameter a of the analytic solution ([T]) so that the 
approximation become accurate enough to study the distribution of the very light modes. 

If the value of C a b c in ([1]) is small outside the range of (|T8|) . the introduction of the cutoff 
will cause small errors. This requires 

aL 2 > 1. (34) 
On the other hand, if C a b c changes gently as a function of momenta, the introduction of 
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spectra 
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Figure 6: The spectra for D — 1, L — 20, a — -p. 

discretization will cause small errors. This requires 

«(P 2 -P' 2 )<1 (35) 
for nearby p, p' . Combining the two conditions, one obtains 

i? ~ a ~ r (36) 

This inequality suggests that the continuum limit L — > oo can safely be taken by simultane- 
ously tuning a in the allowed region. 

The actual process of computing the fluctuation spectra is as follows. The fluctuation 
matrix f[T4"l) is a product of two matrices ^ abc _. As shown in Figure fTTl and fT2l 



def 



ow al 

C=C ' ' dCdE f 



c=c 



contains some loops. The contribution of these loops can be analytically evaluated by Gaussian 
integrations, if the integration variables are the unbounded continuum momenta, but not if 
they are the bounded discrete momenta (TTHT) . In the latter case, the loop integrals must 
be treated numerically. Moreover the existence of the momentum cutoff would also require 
numerical treatment for the internal tree lines. Therefore, to avoid these complexities, I employ 
the unbounded continuous momenta for internal lines, and compute analytically dWabc 

And then the external momenta, i.e. the indices of 



9C de f 

>>w ahc 

c 



c=c 

are assumed to take the bounded 
is diagonalized numerically. This 



9C def 

discretized momenta (fTSj) . and the fluctuation matrix 
procedure would be allowed as an approximation. The analytical expressions of the results of 
the loop integrals in ^ abc _ are given in Appendix O 

Let me first show the result of the numerical computation in D = 1, L = 20, with the 
choice of the parameter a = j2 ■ The spectra are plotted in Figure [6j There one can see 
the clear hierarchy between the modes with the "vanishing" spectra within the resolution of 
the figure and the other "heavy" modes with the spectra > O(0. 1). The number counting of 
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Figure 7: The spectra for D — 2, L = b,a = The horizontal axis is p = a/ (p 1 ) 2 + (p 2 ) 2 - 
The horizontal line is drawn at 0.15. 

the "vanishing" spectra reproduces the result in the preceding section. At p = sector, the 
number is larger by one than the formula (1291) . and agrees with ( l2T?j) at p ^ sectors. This 
implies that there exists a very light mode only at the p = sector. 

Unfortunately, this approximate treatment is not good enough to extract the very light 
mode. In fact, the largest value of the "vanishing" spectra at the p = sector is ~ 3 x 10~ 5 , 
while the second largest is ~ 7 x 10~ 6 , and no clear hierarchy between the zero modes and 
the very light mode can be observed. Note that, in the correct numerical computation in the 
previous section, there existed a clear hierarchy, e.g. ~ 10~ 6 and ~ 10 -15 for D = 1, L = 5. 
Therefore the approximate treatment in this section would be enough only to distinguish 
between the "heavy" modes and the bunch of the zero and the very light modes. However, 
this would be enough to study the distribution of the very light modes by counting the number 
of the bunch of "vanishing" spectra and subtracting the number of the zero modes (I2TJ1) . 

The spectra for D = 2,L = 5, a = jj are plotted in Figure [7J There one finds that the 
"vanishing" spectra form a finite band from the origin of the vertical axis, which would be due 
to the approximate treatment. A hierarchy between the "vanishing" spectra and the other 
"heavy" ones can be clearly observed at p = a/ (p 1 ) 2 + (p 2 ) 2 < 7. One can count the numbers 
of the "vanishing" modes by counting the numbers of modes below a value, say < 0.15. This 
reproduces the observation made in the previous section. Namely, the excess of the number 
in comparison with the formula (1291) is three at p = 0, and one at p ^ 0. 

In the figure, one notices that there exist a series of spectra which leave the "vanishing" 
spectra at p ~ 7 and gradually merges to the "heavy" modes at p ~ 12. This kind of behavior 
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Figure 8: The results for D = 2, L = 5, a = ^ (left), J D = 2,L = 5,a = |f (middle) and 
.D = 2, L = 10, a = jj (right). The last case is studied only for the momentum sectors 
(p\p 2 ) \p l \ < 6, because of the long computational time. 

does not exist in the D = 1 case. These facts suggest that the very light mode begins to mix 
with the "heavy" modes at p > 7. 

The fluctuation spectra for some other parameters are shown in Figure [H The hierarchies 
between the bunch of spectra near the origin and the others can be observed more or less 
clearly in each figure. One can see that, as a becomes smaller, the width of the bunch of 
the "vanishing" spectra tends to become larger, and, as a becomes larger, the "vanishing" 
and "heavy" spectra tend to merge at lower momentum sectors. Therefore, to obtain a clear 
hierarchy between the "vanishing" and the "heavy" modes in a wide range of momentum 
sectors, the parameter a must be chosen appropriately. An appropriate continuum limit 
L — > oo will be obtained with this kind of tuning of a. Counting the number of the spectra in 
each bunch from the origin at low momentum sectors, the same distribution of the very light 
modes as above is obtained. 

In the followings, let us study the higher dimensional cases. In Figure EH the fluctuation 
spectra for D = 3, L = 3, a = are shown. A clear gap between the bunch of spectra from 
the origin and the others can be observed at the momentum sectors p < 2.4. By counting the 
number of the spectra surrounded by the circles in the figure, the excess over the formula 
is obtained as 6 at p = 0, and 3 at p ^ 0. 

In Figure [TU], the fluctuation spectra for D = 4, L = 2, a = jj are shown. The hierarchy 
between the bunch of spectra near the origin and the others is clear for the lowest three 
momentum sectors. By counting the numbers of the spectra surrounded by the circles in the 
figure, the excess over the formula is obtained as 10 at p = and 6 at p ^ 0. 
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Figure 9: The fluctuation spectra for D = 3, L = 3, a = ^jj-. The numbers of the spectra near 
the origin surrounded by the circles are counted and compared with the formula fl29l) . 
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Figure 10: The fluctuation spectra for D = 4, L = 2, a = With this small value of L, 
the hierarchy between the bunch of spectra near the origin (the spectra surrounded by circles) 
and the others is clear only for p — 0, 1, \f2. 
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D 


flight modes at p = 


# light modes at p ^ 


1 


1 





2 


3 


1 


3 


6 


3 


4 


10 


6 



Table 1: The distribution of the very light modes obtained from the numerical analysis. 

6 A physical interpretation of the very light modes 

In the preceding two sections, it is observed that the distribution of the very light modes has 
a characteristic feature depending on the dimensions, which is summarized in Table [TJ A 
general expectation in lattice-like approaches is that, in the continuum limit L — > oo, the low- 
energy effective physics is governed by such light modes, while the "heavy" spectra > O(0. 1) 
appear only at a "Planck" energy, and the zero modes are just the unphysical gauge modes. 
In this section, it will be pointed out that the distribution of the very light modes can easily 
be explained by the general relativity. 

The fundamental variable of the general relativity is a metric tensor. Since the general 
relativity has the gauge symmetry of the general coordinate transformation, the real degrees 
of freedom are the metric tensor modulo the general coordinate transformation. These real 
degrees of freedom should be compared with the numerical result. As mentioned previously, 
the background space corresponding to the solution with ([18]) is a fuzzy D-dimensional 
flat torus. The infinitesimally small general coordinate transformation around a flat space 
background is given by 

&9nv = d^v u + d u v^ (37) 
where is the local translational transformation vector. In the momentum basis, this is 

Sg^uip) = W»v v {p) + ipvV^p). (38) 

Therefore, in the p = sector, the gauge transformation does not generate any transforma- 
tions, while, in the p / sectors, it generates a D dimensional transformation. Thus the 
degrees of the freedom of the metric tensor modulo the general coordinate transformation are 
given by 

f EiE+R f or „ = o 
#Wg.ct. = ( ^ for ^ ; (39) 

This formula exactly agrees with Table [U 
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7 Discussions and future prospects 



In this paper, numerical analysis of the fluctuation spectra around a Gaussian classical solution 
of a tensor model was carried out. It was observed that the distribution of the very light modes, 
which are expected to be the low-energy effective physical modes in the continuum limit, 
agrees exactly with what is expected from the general relativity in one to four dimensions. 
This result suggests that the low-energy effective field theory of the tensor model around the 
classical solution can be described in a similar manner as the general relativity. 

The exact agreement in various dimensions would be a good evidence of the link between 
the tensor model and the general relativity, but it is clear that more evidences are needed 
to definitely prove it. A direction would be to study the properties of the very light modes, 
and compare them with the general relativity. For this purpose, the very light modes must be 
distinguished from the other modes, especially from the zero modes. This can be done with the 
correct numerical method employed in Section HJ but not with the efficient but approximate 
method in Section O since, in the latter method, the very light modes and the zero modes 
form a bunch. Therefore, with the present computational facility mentioned in Section HI the 
properties of the very light modes can be studied only in one and two dimensions for small L. 

To study the properties of the very light modes in larger D, L, it would be necessary to 
develop a systematic method to reduce the number of the degrees of freedom of the tensor 
model without affecting the properties of the very light modes. This would also be necessary to 
study the quantum properties of the tensor model. One candidate for such a method would be 
a kind of renormalization procedure for the tensor model. A key issue would be the treatment 
of the zero modes, which can easily mix with the very light modes in any rough approximate 
method like that in Section [5j It would be important to keep the symmetry of the tensor 
model, so that the zero modes be distinguishable from the very light modes. 

Besides the huge number of the degrees of freedom of the dynamical variables, what makes 
tensor models quite complicated is the lack of a guiding principle to chose among the infinite 
possibilities of actions. It would be necessary to single out the elements essential to induce 
the general relativity at low energy, and get a hind for a principle. It would be also important 
to investigate the universality classes of tensor models, to see whether the result of this paper 
is just a special case or universal in tensor models. 

Finally it would be amusing to see that some constructive approaches to string/M theory 
[531 151] have dynamical variables (M^) a b , which is similar to a rank-three tensor. A tensor 
model is discussed also as a toy model of open membrane theory [20]. Thus, as explained 
in Section [TJ tensor models appear in various approaches to quantum gravity, i.e. dynamical 
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triangulation, loop quantum gravity, fuzzy spaces and string/M theory. It would be highly 
interesting, if tensor models connect all of these apparently different approaches in future. 
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A The computation of the partial derivatives 

In this appendix, the explicit formula of the partial derivatives in (|14|) are given. 
The partial derivative of the first term of W a b c in ([6]) is obvious, 

Here the permutation symmetry of the indices of the tensor C a b c will be cared by the summation 
over the permutations a, a' in (|14p . The same kind of simplification of expressions will also 
be used in the formulas below. 

The partial derivative of the second term of W a b c consists of two parts. One is the derivative 
a ° 6 . This can be done through the well known formula, 

J^T- = ~ K aa' 7^ K Vb ■ ( 41 ) 

The definition of K a b is given by neglecting the bar ~ in ([3l). From this definition, the derivative 
7p— can be obtained as the diagram in Figure [TTJ 

The other part is the partial derivative J^y- The definition of I a b c is given by neglecting 
the bar ~ in dSJ). The result of the partial derivative is shown in Figure [T2l In this figure, to 
avoid a large figure, it is omitted to average over all the permutations of the open ends in the 
left, which is needed to take into account the permutation symmetry of the indices of I a b c . 
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Figure 11: The partial derivative ar ab . The open ends in the left are associated with the 
indices of K ab} and those of C c a e are in the right. The permutation symmetry of the indices 
of C a bc will be cared by the summation of a, a' in (THj) . 




indices of I a bc, and those of Cdef are in the right. The permutation symmetry of the indices of 
lohc must be cared by averaging over all the permutations of the open ends in the left, while 
the permutation symmetry of the indices of C a b c are cared by the summation of a, a' in (fT4l) . 
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-3 3 


-3 1 2 


-2 2 


1 .84462617663699010E-001 


2.40366883259875430E-001 


3.61070116973867590E-001 


-2 1 1 


-1 1 





4.17490143455173490E-001 


5.58142487944932600E-001 


6.53570115178405460E-001 



Table 2: The solution C piP2P3 for D = 1, L = 3. The momenta P1P2P3 are listed in the upper 
colums, and the corresponding values of C° p3 are in the lower. 

B The raw numerical results 

B.l The numerical results for D = 1, L = 3 

Because of the permutation symmetries of the indices of C P1P2P3 and the ansatz, there are only 
6 independent variables in the optimization process. The minimum of ( 13T1) reached by the 
Nelder-Mead method is f{C) ~2x 10 -25 . The values of the six C p at the minimum are 
listed in Table EJ 

In Table [3l the fluctuation spectra obtained by diagonalizing (fT4l) are shown for each 
momentum sector. The spectra contains very tiny vales < 10 -14 . These can be identified with 
the zero modes from the symmetry breaking SO(n, R)/SO(2), since the numbers agree with 
the formula (1291) . This agreement is also a good check for the correctness of the computational 
codes, because any slight mistakes in the codes will break this agreement. 

There exists a non-zero but a very small spectrum ~ 2 x 10~ 3 only at p = sector. The 
other spectra are > 0(0.1). 

B.2 The numerical results for D = 1, L = 5 

The minimum value reached is f(C) ~ 3 x 10~ 25 . The solution is shown in Table HI The 
fluctuation spectra at the p = sector are shown in Table El There the tiny spectra < 10~ 15 
can be identified with the zero modes, the number of which agrees with the formula (129]) . 
There exists a very light mode ~ 10~ 6 , and the others are > O(l)- The fluctuation spectra at 
the other momentum sectors are also obtained, but the values are not listed here. The spectra 
are plotted in Figure [3] instead. At the p 7^ sectors, the fluctuation spectra can be divided 
into two parts. One part consists of the tiny spectra < 10~ 13 , the numbers of which agree 
with (I2"9"j) . The other consists of the " heavy" modes with the spectra > O(0.1). 
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V 


spectra 


#zero 




-2.074144E-015 1.296988E-015 2.010206E-003 1.187204E+000 







2.226274E+000 6.450744E+000 2.101532E+001 1.343863E+003 


2 




-8.779625E-014 -1.021941E-014 -1.835845E-015 1.208359E+000 




1 


1.623874E+000 3.193161E+000 2.674381E+001 1.130220E+003 


3 




-5.557866E-015 5.578683E-015 1.018426E+000 1.622450E+000 




2 


3.177940E+000 3.086634E+001 6.450256E+002 


2 




-a -a a r\ >~~7 r\ v t — \ r\ a a o r~\ c~\ A A v a t — \ r\ a /~* /~\ r\ a r\ t~\ A t — \ r\ r\A -a n v a >~~j r\ r\~\ — \ > r\r\c\ 

-1.110795E-014 -3.234454E-016 9.040864E-001 1.751722E+000 




3 


2.161799E+000 2.207987E+001 2.123835E+002 


2 




-5.834977E-015 5.523345E-001 1.344oUoE+UUU 4.7o3487E+000 




4 


A A^/IQT t^T? I AAI 

4.Ub4ol0il/+UUl 


1 


5 


1.792788E-016 3.732453E-001 1.315828E+000 1.723215E+001 


1 


6 


6.317254E-001 9.912976E-001 7.28901 1E+000 





7 


4.695151E-001 1.014654E+000 





8 


9.269276E-001 





9 


9.965088E-001 






Table 3: The values of the fluctuation spectra obtained by diagonalizing (Ti"4"|) for D = 1, L — 3. 
The momentum of the modes, the fluctuation spectra, and the number of the zero modes from 
the formula (1291) are listed for each momentum sector. The negative p sectors are omitted as 
obvious, since the ansatz assumes the reflection symmetry p — > —p. 



-5 5 


-5 1 4 


-5 2 3 


5.30666376133652470E-002 


7.72763400282020090E-002 


9.34623815478460020E-002 


-4 4 


-4 1 3 


-4 2 2 


1.24046348049658310E-001 


1.65453909450288150E-001 


1 .82235408294834920E-001 


-3 3 


-3 1 2 


-2 2 


2.43404073757127580E-001 


2.95609691797031670E-001 


3.95946657504184290E-001 


-2 1 1 


-1 1 





4.36599185858296160E-001 


5.30865955626977360E-001 


5.85474131045624180E-001 



Table 4: The table of the solution C® for D = 1, L = 5. 



V 


spectra 


^zero 





-2.796163E-015 -8.151607E-016 2.246353E-016 1.433114E-015 
1.073950E-006 9.152090E-001 9.941203E-001 1.015528E+000 
1.037451E+000 1.116800E+000 1.242951E+000 1.383553E+000 
1.934481E+000 4.332078E+000 5.159299E+000 8.182407E+000 
2.407250E+001 1.363542E+003 


4 



Table 5: The fluctuation spectra in the p = sector for D = 1, L = 5, and the number of the 
zero modes from (129]) . 
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(-2,-2) (0,0) (2,2) 


(-2,-2) (0,1) (2,1) 


(-2,-2) (0,2) (2,0) 


1 .33090366933366510E-001 


1.54486879640217570E-001 


1 .33090366933252490E-001 


(-2,-2) (1,1) (1,1) 


(-2,-1) (0,0) (2,1) 


(-2,-1) (0,1) (2,0) 


1.79323241274937400E-001 


2.04924589321128730E-001 


2.04924589320929830E-001 


(-2,-1) (0,2) (2,-1) 

\ 7 / \ 7 / \ 7 J 


(-2,-1) (1,-1) (1,2) 


(-2,-1) (1,0) (1,1) 


1.54486879639913640E-001 


1.79323241274623960E-001 


2.37869660255753430E-001 


(-2,0) (0,0) (2,0) 


(-2,0) (1,-1) (1,1) 


(-2,0) (1,0) (1,0) 


2.54117050163986060E-001 


2.37869660255717370E-001 


2.94970635725226470E-001 


(-1,-1) (0,0) (1,1) 


(-1,-1) (0,1) (1,0) 


(-1,0) (0,0) (1,0) 


3.15530629872480820E-001 


3.15530629872307290E-001 


3.91274239784655140E-001 


(0,0) (0,0) (0,0) 






4.85200219013368070E-001 







Table 6: The solution C? i 2u i 2m i 2n for D = 2, L = 2. 
B.3 The numerical results for D = 1, L = 10 

The number of independent variables is 36, and the minimizing process takes a longer time. 
A one-day run was carried out, but a minimum vanishing within the computational accuracy 
was not reached. The minimizing process was terminated at f{C) ~ 10~ 14 . The fluctuation 
spectra are plotted in Figure [31 The qualitative features are the same as L = 3, 5. A clear 
agreement between the number of the tiny spectra and (fT4"|) is obtained. There exists again a 
very light mode only at the p = sector, which has a spectrum ~ 3 x 10~ 7 . 

B.4 The numerical results for D = 2, L = 2 

The minimum value reached is f(C) ~ 10~ 25 . The values of C^ bc at the minimum are shown in 
Table [6j The fluctuation spectra at all the momentum sectors are computed by diagonalizing 
frT4|) . To avoid a huge table, only the lower parts of the spectra in the p = (0, 0) and p = (1, 0) 
sectors are shown in Tabled There the number of modes with values < 10 -14 agree with the 
formula (129]) . and can be identified with the zero modes. There seem to exist a number of very 
light modes. In the p = (0, 0) sector, the three modes with ~ 8 x 10~ 9 , 3 x 10~ 3 , 3 x 10 -2 may 
be identified with the very light modes. In the p = (1,0) sector, the mode with ~4x 10~ 2 
may be identified with a very light mode. The identification of the very light modes is not 
so clear in this case, but will be improved much more in L = 3. The spectra in all the other 
momentum sectors are computed as well, but are not shown here. There the modes with 
spectra < 10 -13 can clearly be identified with the zero modes, and their numbers agree with 



24 



p 


spectra 


#zero 


(0,0) 


-3.422193E-015 -2.176675E-015 -8.331255E-016 -2.168759E-016 
5.312991E-017 1.245100E-016 6.100288E-016 1.548574E-015 
3.822409E-015 6.586865E-015 8.393766E-009 2.734072E-003 
2.650856E-002 1.375314E-001 5.215319E-001 5.803728E-001 
5.930927E-001 7.295304E-001 • • • 


10 


(1,0) 


-2.247079E-014 -6.834805E-015 -3.502970E-015 -1.232192E-015 
-6.677623E-016 6.984656E-017 1.852255E-016 3.004076E-015 
5.622760E-015 1.195054E-014 4.019257E-002 2.451403E-001 
2.749515E-001 5.590763E-001 5.978878E-001 • • ■ 
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Table 7: The lower part of the fluctuation spectra at the p — (0, 0), (1, 0) sectors for D = 2, 
L = 2, and the number of the zero modes fT2"9~l) . 

the formula (1291) . There the very light modes and the modes with spectra > O(0.1) cannot be 
well distinguished. 



B.5 The numerical results for D = 2, L = 3 

The number of independent variables is 42. The minimum value reached is f(C) ~ 10~ 22 . To 
avoid a large table, the solution is not explicitly shown in this case. The fluctuation spectra 
are computed in all the momentum sectors. Some parts of the spectra in some low momentum 
sectors are shown in Table [SJ In the p = (0, 0) sector, the formula (1291) implies that the modes 
with spectra < 10 -14 should be regarded as the zero modes. Then the very light modes can 
be identified with those with the spectra 1.3 x 10 -12 , 3.0 x 10~ 4 and 1.1 x 10~ 3 . In the non- 
vanishing momentum sectors, the number of the modes with spectra < 10~ 13 agrees with the 
formula ([29]). 



C Analytical expressions of the loop integrals in 



dW, 



abc 



dC def 



dW ab 



c=c 



In this appendix, the analytical results of the loop integrals necessary for evaluating , . 
are given. The diagrams in Figure [TT] and [TJ] show that the kinds of loop integrals can be 



summarized as in Figure | 

The explicit analytical expressions for these diagrams are given in the followings. The 
expressions must be multiplied by the obvious delta functions of momentum conservation, 
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(0,0) 


-7.023228E-015 -3.339147E-015 -1.354682E-015 -1.152946E-015 
-1.022040E-015 -8.653429E-016 -5.793971E-016 -4.444370E-016 
-2.825529E-016 -1.526719E-016 -9.835320E-017 7.877653E-017 
2.207186E-016 3.057724E-016 4.016902E-016 5.873527E-016 
8.141446E-016 1.301374E-015 1.455441E-015 2.779198E-015 
5.143835E-015 1.357294E-014 1.312193E-012 3.048268E-004 
1.069215E-003 1.385427E-001 3.174581E-001 ••• 


(1,0) 


• • • 3.531105E-014 2.117237E-003 2.150369E-001 • • • 


(1,1) 


• • • 2.769432E-014 9.569329E-004 1.232058E-001 • • • 


(2,0) 


• • • 6.895682E-015 6.938484E-003 2.403594E-001 • • • 



Table 8: A part of the fluctuation spectra for D = 2, L = 3. For the p = (0, 0) sector, all the 
lower spectra are shown. For p = (1, 0), (1, 1), (2, 0) sectors, only the very light modes and the 
nearby spectra are shown. 




26 



which are omitted below for simplification, e.g. S (pi + p 2 + P3) for Gl, etc. 



Gl 
G2 



G3 



G4 



G5 



A 



7T \ T 



6a J 



cxp 

ID 



A 



36 



A 



^ (-3pi + 3p 2 - 7p 3 + 7p 4 ) 2 + 3 (p a + p 4 ) 2 + 3 (pi + p 3 )' 

r. 



41 



7T 



a 



8D 



7T 



7T 



18a/ V26a 



exp [-a (3pi + 3^2 + 5p 3 + lp\ + p 2 , 



13 



-3pi - 5p 3 + 5p 4 + 2p 5 ) 2 + 4(p 4 + p 5 ) 2 + 6(pi + p 3 )' 



44 



V4v / 2^ 



8.D 



7T 



18a 



1? 



7T 



4 ■ 7 • 13a 5 



cxp 



(749p 2 + 82 \p\ + 42 lp 2 + 50p 3 p 4 

91 



(42) 



(43) 



(44) 



+ 177p 2 4 + 132p 3 p 5 + 10p 4 p 5 + I77p\ + 396p 3 p 6 + 30p 4 p 6 + 152p 5 p 6 + 319p^+ 
14 Pl (47p 2 + 15p 3 + 10p 4 + 3p 5 + 9p 6 ) + 2p 2 (150p 3 + 61p 4 + 30p 5 + 90p 6 ))] . 



(45) 



.4 



11 



7T 



Ay/2. 



a 



8D 



71 



18a 5 



D 



8-47 



exp 



a 



(l499p 2 + lOllp 2 , + 1283p 2 + 370p 3 p 4 



+803p^ + 222p 3 p 5 + 738p 4 p 5 + 635^ + 74p 3 p 6 + 246p 4 p 6 + 298p 5 p 6 + 363p^+ 

10p 2 (41p 3 + ll(5p 4 + 3p 5 + Pe)) + 2pa (313p 2 + 443p 3 + 73(5p 4 + 3p 5 + p 6 )))] . (46) 
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